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Using first-principles approaches, this study investigated the efficiency of energy conversion in 
nanojunctions, described by the thermoelectric figure of merit ZT. We obtained the qualitative and 
quantitative descriptions for the dependence of ZT on temperatures and lengths. A characteristic 
temperature: To = vWtW was observed. When T < To, ZT oc When T > To, ZT tends to a 
saturation value. The dependence of ZT on the wire length for the metallic atomic chains is opposite 
to that for the insulating molecules: for aluminum atomic (conducting) wires, the saturation value 
of ZT increases as the length increases; while for alkanethiol (insulating) chains, the saturation 
value of ZT decreases as the length increases. ZT can also be enhanced by choosing low-elasticity 
bridging materials or creating poor thermal contacts in nanojunctions. 



There has been renewed interest in the study of ther- 
moelectricity motivated by its possible application in 
energy-conversion devices at the nanoscale level [l|, 0, 
S, i, i, @, 0, B B E ini- Recent experiments on the 
Seebeck coefficient, which is insensitive to the number 
of molecules in the junction, shed light on the possibil- 
ity of the implementation of thermoelectric devices at 
the atomic level Q. Nanoscale energy-conversion de- 
vices can convert waste heat energy into useful electric 
power and stabilize miniature electronic devices by re- 
ducing the temperature. The Seebeck coefficient, which 
is related not only to the magnitude but also to the slope 
of density of states, can provide more information than 
current- voltage characteristics P, d, H, Q • Measurement 
of the Seebeck coefficient has been applied to explore the 
effect of chemical structure on the electronic structure of 
molecular junctions [1]. The gate field has been theoreti- 
cally proposed as a means of modulating the conduction 
mechanism between p-type [the Fermi energy is closer to 
the highest occupied molecular orbital (HOMO)] and n- 
typc [the Fermi energy is closer to the lowest unoccupied 
molecular orbital (LUMO)] via the sign of the Seebeck 
coefficient d, |Tl[ . Although much research has been de- 
voted to the study of Seebeck coefficient, little is known 
about the efficiency of energy conversion in nanojunc- 
tions [1]. The objective of this research was to provide 
greater insight into this subject. 

Molecular tunneling junctions consist of source-drain 
electrodes as independent electron and heat reservoirs 
with distinct temperatures [T2^(fl)] and chemical poten- 
tials [^JiL(R)\■ The efficiency of energy conversion depends 
on several factors: the electrical conductance (cr), the 
Seebeck coefficient (S*), the electron thermal conductance 
(kez) and the phonon thermal conductance {Kph)- The 
efficiency can be described by the dimensionless thermo- 
electric figure of merit [^: 



ZT : 



-T, 



(1) 



where T = {Tl+Tr) /2 is the average temperature in the 
source-drain electrodes. The ideal thermoelectric molec- 
ular junction would have a large S, a large u and a small 



combined thermal conductance (keZ + Uph)- Thermoelec- 
tric materials with a large a are usually accompanied by 
a large k^i , which makes the enhancement of the thermo- 
electric figure of merit a challenging task. 

In this Letter, we reported first-principles calculations 
of the thermoelectric figure of merit in nanojimctions. It 
aimed to obtain a qualitative and quantitative descrip- 
tions of ZT for temperatures and lengths of the nano- 
junctions. The self-consistent density functional theory 
(DFT) was performed together with the derivation of an 
analytical expression for ZT to investigate its dependence 
on the temperatures and lengths of nanojunctions. As an 
example, this study investigated ZT for the aluminum 
atomic (conducting) wires and the alkanethiol (insulat- 
ing) molecules in a nanojunction in the linear response 
regime. It was found ZT oc T^ at low temperatures, 
while ZT tended to a saturation value at high tempera- 
tures. The dependence of ZT on the wire lengths for the 
metallic atomic chains was opposite to that for the in- 
sulating molecules: longer conducting wires and shorter 
insulating molecules had better efficiency of energy con- 
version. The results of this study may be of interest 
to experimentalists attempting to develop thermoelectric 
nanoscale devices. 

First, we started by a brief introduction of the DFT 
calculations for a molecule sandwiched between two bulk 
electrodes with external source-drain bias. The effec- 
tive single-particle wave functions of the whole system 
were calculated in scattering approach by solving the 
Lippmann-Schwinger equation with exchange and corre- 
lation energy included within the local density approxi- 
mation itcratively until the self-consistency was obtained. 
The effective single-particle wave function *^^"^(r,K|[) 
represents the electron incidents from the left (right) elec- 
trode with the energy E and component of the momen- 
tum K|j parallel to the electrode surface [H, El, fl^ - 
These left- and right-moving wave functions, weight- 
ing with the Fermi-Dirac distribution function accord- 
ing to their energies, were applied to calculate the elec- 
tric current / and the thermal current conveyed by the 
transport electrons Jq , via the following expressions: 
I=fJdE [f§{^iR, Tn)T^{E) - fk{t^L,TL)r^{E)] and 
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= f / dE [{E - f§r"iE) - {E - /|r^(i?)] , 



where t^^^\E) is the 
the electron with energy 
(right) electrode. The 
be computed using the 
self-consistently in DFT 
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transmission function of 
E incident from the left 
transmission function can 
wave functions obtained 
calculations according to 
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and dR 



represents an element of the electrode surface. It is as- 
sumed that the left and right electrodes served as inde- 
pendent electron and phonon reservoirs with the electron 
population described by the Fermi-Dirac distribution 



function, /^'^^ = 1/ (cxp {{E - p.L(R)) IkBTL(R)) + l), 
where ks is the Boltzmann constant, and ^.l{r) a-nd 
T]^(R) are the chemical potential and the temperature 
in the left (right) electrode, respectively. The external 
source-drain bias is defined by: Vb = [i^r. — A'l)/^- 

Then, we briefly described the method used to cal- 
culate the electrical conductance, the Seebeck coefficient 
and the thermal conductance conveyed by electron trans- 
port. We considered the extra electric and thermal cur- 
rent induced by an additional infinitesimal temperature 
(AT) and voltage (Al/) symmetrically distributed across 
the junction: 



J 



AT AT eAV eAV 

A/ = I{^il.Tl + —;fiR,TR-—) + I{fiL + -^,TL;fiR - -^,Tr) - 2/(^i, TL;fiR, Tr), (2) 



and 



AT AT 



^^^,Tl;hr - ^^,Tr) - 2Jq{^.l,Tl]^.r,Tr), (3) 



respectively. After expanding the Fermi-Dirac distribu- 
tion function to the first order in AT and Ay, we ob- 
tained the Seebeck coefficient (defined by S" = AF/AT) 
by letting A/ = and the electron thermal conductance 
(defined by k^i = AJ^/AT): 



S 



^ Tl ^ Tr 



Kel 
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(4) 
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where if,^'^' 



jdE (E-fiLiR)) 



dE 



-t{E), and 



t{E) = T (E) = T (E), a direct consequence of the 
time-reversal symmetry. In addition, the differential con- 
ductance, typically insensitive to temperature in cases 
where direct tunneling is the major transport mechanism, 
may be expressed as: 



E 

i=L,R 



ksTi 



■^'^K{E)dE. 



(6) 



So far, the physical quantities that have been discussed 
have been related to the propagation of electrons. How- 
ever, in most cases, the thermal current is dominated 
by the contribution from phonon transport. In the ab- 
sence of the phonon thermal conductance, the research 
on ZT is incomplete. To consider the phonon contribu- 
tion to ZT, it is assumed that the nanojunction is a weak 
elastic link, with a given stiffness that may be evaluated 
from total energy calculations, attached to the electrodes 



modeled as phonon reservoirs. We estimate the contribu- 
tion of the thermal current from phonon scattering (Jn^ ) , 



following the approach of Patthon and Geller [15|. Af- 
ter expanding the Bose-Einstein distribution function to 
the first order of AT in the expression of phonon ther- 
mal current, the phonon thermal conductancc(defined by 
kph = AJq'/AT) is obtained: 



hkn 



f dEE^NLiE)NRiE) E ■ 

i=L.R 



n,{E){l + n,{E)) 



(7) 

where nL(R) = l/(e^/^B^^(«) - 1) and NL(R)iE) ~ CE 
is the Bose-Einstein distribution function and the spec- 
tral density of phonon states in the left (right) electrode, 
respectively. The stiffness of the bridging nano-structure 
is: K = YA/l, where Y is the Young's modulus and A 
(l) is its cross-section (length). 

Finally, ZT could be calculated by applying Eqs. (|4]) 
to (O. The Seebeck coefficient and the electron (phonon) 
thermal conductance can be characterized by the power 
law expansions: S w aT, k^i ~ (3 [T + rjT^] w /3T 
and Kph = j{l)T^ in the common range of tempera- 
tures (Ti « Tr = T) and in the linear response regime 

{p-L « /i_R = m): where a = -Tr^^B^gB^/ (3er(//)); 
a = 27r2fc|r(/.)/(3/i); ^ = {i^ksdriy.) j dEf j {Zr^p)'') 
and 7(0 = %-K^k%C'^A^Y'^ lixm"^). Consequently, the 
thermoelectric figure of merit in the nanojunctions has a 
simple form. 



ZT: 



a' 



■aT^ 



/3T + 7(0T3 



(8) 
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FIG. 1: Aluminum atomic junctions at Vs = 0.01 V: (a) 
Schematic of 3-Al atomic chain and its Energy diagram. The 
Al-Al bond distance was about 6.3 a.u.; (b) Electrical con- 
ductances a vs T; (c) Electron thermal conductances k^i vs 
T; (d) Seebeck coefficients S vs T; (e) Phonon thermal con- 
ductances vs T (F = 1.2 X 10^^ dyne/cm^); (f) Phonon 
thermal conductances vs T (F = 5.0 x 10^ dyne/cm^); 
(g) Log{ZT) vs T (for F = 1.2 x lO" dyne/cm^); (h) ZT and 
Log{ZT) (Inset) vs T (for F = 5.0 x 10^ dyne/cm^). 



which is vahd in small bias and low temperature regimes. 

The properties of the thermoelectric figure of merit 
now can be discussed using Eq. ([S]). There was a char- 
acteristic temperature, Tq = for ZT in the 
nanojunctions. When T ^ Tq, the thermal current was 
dominated by the contribution from the electron trans- 
port {kph ^ fcei), which led to ZT increasing as the 
temperature increased: ZT w aS^T/kei ~ [a'^a/0\ T^. 
Similarly, when T ^ Tq, the thermal current was dom- 
inated by the contribution from the phonon transport 
{kph ^ fcez), which led to a saturation of ZT at a 
constant value related to the length of the junction: 
ZT « aS'^T/kp « aV/7(0. To increase ZT it was 
first necessary to reduce kph by choosing low-elasticity 
bridging wires or creating poor thermal contacts in the 
nanojunctions, such that ZT « aS'^T /kei- It is worth 
noting that a and k^i roughly canceled each other out in 
the contribution of ZT because both were proportional 
to T(/Li). It then followed that ZT oc S'^T and, thus, that 
the material with a large Seebeck coefficient was of key 
importance to increasing ZT . The characteristic mark 
of such a material in the nanojunctions is a sharp peak 
around the Fermi levels in the DOS [llj, and the Seebeck 
coefficient may be optimized by applying the gate field 
[3, UH . In addition, it was noted that a and a depend on 
the length of the junction in a way related to the mate- 



FIG. 2: Alkanethiol junctions at Vb = 0.01 V: (a) Electric 
conductance a vs T; (b) Electron thermal conductances k^i 
vs T; (c) Seebeck coefficients S vs T; (d) Phonon thermal 
conductance vs T (F ~ 2.3 x 10^^ dyne/cm^); (e) ZT 
vs T (F ~ 2.3 x 10^^ dyne/cm^); (f) ZT vs T (for F = 
dyne/cm^). 



rial properties of bridging wires, which is refiected in the 
distinguished features of ZT on the length dependence. 
This point was explained using two catalogs of nanojunc- 
tions: the aluminum atomic (conducting) wires and the 
alkanethiol (insulating) chains, as discussed below. 

Aluminum atomic wire is ideal for studying charge 
transport at the atom-scale [see Fig . |T[a) for a schematic 
of the aluminum junction] |l6l . 17. Tisl. [l^. As shown 
in Fig. mb), the conductance was relatively insensitive 
to the chain length (typically around 1 Go = 2e^//i « 
77 ^S) apart from the possible 4-atom periodicity due 
to a filling factor of 1/4 in the tt orbitals [i^]- As shown 
in Fig. [ijc), the magnitude of electron thermal conduc- 
tance was linear in temperatures, Kg; ~ [3T . At a fixed 
temperature, the dependence of the magnitude of Kg; on 
the number of Al atoms was the same as that of cr, ow- 
ing to the fact that both a and n^i were proportional 
to r(/i). As shown in Fig. [IJd), the magnitude of the 
Seebeck coefficient was linear in temperature, S k, oT, 
with the negative sign showing that the carrier was n- 
type. At a fixed temperature, it was observed that the 
magnitude of the Seebeck coefficient increased consid- 
erably as the number of Al atoms increased. The in- 
crease of the Seebeck coefficient was due to the increase 
of the slope in the DOS at the Fermi level. These fea- 
tures may be related to the fact that the Fermi level 
was close to the LUMO in the Al wires. Fig. [Ije) shows 
the phonon thermal conductance: Kph = j{l)T^, for the 
Young modulus using Y ~ 1.2 x 10^^ dyne/cm^ from 
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the total energy calculations 2l|. As seen, Kph ^ KeZ 
was due to the large Young modulus. Fig.[ljg) shows the 
thermoelectric figure of merit with Uph calculated using 
Y = 1.2 X 10^^ dyne/cm^ from the total energy calcu- 
lations. The increase in the number of Al atoms sharply 
increased the saturation value of ZT because of the sharp 
increase in the Seebeck coefficient by the number of Al 
atoms according to ZT oc S^. The thermoelectric figure 
of merit reached the saturation value, ZT —> a'^a/j{l) 
when T ^ Tq. Since the mechanical elasticity of the Al 
wires could be delicate to the detailed geometry in the 
contact region which was unknown in the real experi- 
ment, as Fig.[T]^f) shows the phonon thermal conductance 
Kph for another possible value, Y — 5.0 x 10^ dyne/cm^ 
from controlled tensile experiments on nanoscale Al films 
(2^ . As shown in Fig. [TJh), it is worth noting that ZT 
could be strongly enhanced by a smaller Kph- In such 
cases, the thermal current conveyed by electron trans- 
port dominated so that ZT [^o?a/ (3) T^ and ZT was 
strongly enhanced. 

Alkancthiols [CH3(CH2)„-iSH, denoted as C„] are a 
good exam ple of reproducible junctions that can be fab- 
ricated 13, IS] ■ In contrast to the conductor behavior of 
aluminum wires, alkanethiol chains are insulators. It has 
been established that non-resonant tunneling is the main 
conduction mechanism in alkanethiol junctions. Conse- 
quently, the conductance is small and decreases expo- 
nentially with the length of wire, as a = ao exp {—^l) 
where I is the length of alkanethiol chain and ^ « 0.78 
A-i m, M, [13, [H, [li], as shown in Fig. M^a). By ex- 
ploiting the periodicity in the {CH2)2 group of the alka- 
nethiol chains, the wave functions of the C„ junctions 
were calculated by a simple scaling argument, which led 
to exponential scaling in the transmission function t{E). 
As shown in Fig. [^b), the magnitude of electron ther- 
mal conductance was linear in temperatures, n^i ~ f3T. 
At a fixed temperature, the magnitude of k^i decreased 
exponentially with n, the number of carbon atoms in 
Cn, owing to the scaling behavior of t{E). As shown in 
Fig. IHc), the magnitude of the Seebeck coefficient was 
linear in temperature as 5' « — Tr^fc^ / (3eT(/x))r, 
and its dependence on the number of carbon atoms was 
canceled due to the same scaling factor exp {—£,1) for 
both t(/z) and ^g^. As shown in Fig. HJc), thermal 
conductance increased as the temperature increased as 
f^ph ~ l{l)T^ for the Yo ung modulus calculated with 
total energy calculations [2l|. At a fixed temperature, 
Kpii decreased as due to oc (see Fig. [2ljd)). 
Due to the small transmission probability for the insulat- 
ing alkanethiol chains, the electron thermal conductance 
(note: kei cc cr) was much suppressed so that Kei ^ kph, 
as shown in Fig. [5]^b) and (d). Consequently, the charac- 
teristic temperature Tq was low in the alkanethiol chains, 
and the T^ regime for ZT was significantly suppressed. 
As shown in Fig. Elje). ZT decreased as the number of 
carbon atoms increased for T ^ Tq, due to the satu- 
ration value of ZT « a^a/^{l) cc Pexp(— ^/). Never- 



theless, there was enough experimental evidence to show 
that the junctions had poor thermal contacts for certain 
samples [30|. These samples quickly frustrated at much 
smaller biases. The frustration of these samples could 
have been due to poor heat dissipation by the thermal 
current via phonon transport when the local heating was 
triggered by an external bias larger than the threshold 
value IJI, ISlJ. In such cases, the bridging nano-structure 
effectively has a very small Young modulus. In the limit 
of extremely poor thermal contacts (effectively, kph = 0, 
cr and kei canceled the length dependence), which led to 
ZT cc a^aoT^/P and impfied that ZT was independent 
of the wire length as shown in Fig. El^g). 

In conclusion, self-consistent DFT calculations were 
performed to study the efficiency of energy conversion 
in nano junctions. There was a characteristic tempera- 
ture of To for ZT: when T < Tq, ZT (a'^a/P) T^; 
when T > Tq, ZT a^a/^{l). Of key importance to 
increasing the efficiency of energy conversion was using 
materials with a large Seebeck coefficient. Such materials 
were usually characterized by a sharp peak around the 
Fermi levels in the DOS. Efficiency could be further opti- 
mized by applying the gate field, choosing low-elasticity 
bridging materials or creating poor thermal contacts in 
nanojunctions. The relation between ZT and the wire 
lengths depended on the material properties: for alu- 
minum atomic (conducting) wires, the saturation value 
of ZT increased as the length increased; while for the 
alkanethiol (insulating) chains, the saturation value of 
ZT decreased as the length increased. The conclusions 
of this study may be beneficial to research attempting 
to increase the efficiency of energy conversion in nano 
thermoelectric devices. 
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